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Abstract 

Aims. We wish to study the spectral dependence of the radio emission from cosmic-ray air showers around 100 PeV (10 17 eV). 
Methods. We observe short radio pulses in a broad frequency band with the dipole-interferometer LOPES (LOFAR Prototype 
Station), which is triggered by a particle detector array named Karlsruhe Shower Core and Array Detector (KASCADE). LOFAR 
is the Low Frequency Array. For this analysis, 23 strong air shower events are selected using parameters from KASCADE. The 
radio data are digitally beam-formed before the spectra are determined by sub-band filtering and fast Fourier transformation. 
Results. The resulting electric field spectra fall off to higher frequencies. An average electric field spectrum is fitted with an 
exponential E v = K ■ exp(^/MHz//3) and (3 = —0.017 ± 0.004, or alternatively, with a power law e v = K ■ v a and a spectral 
index of a = — 1 ± 0.2. The spectral slope obtained is not consistent within uncertainties and it is slightly steeper than the slope 
obtained from Monte Carlo simulations based on air showers simulated with CORSIKA (Cosmic Ray Simulations for KASCADE). 
For the analyzed sample of LOPES events, we do not find any significant dependence of the spectral slope on the electric field 
amplitude, the azimuth angle, the zenith angle, the curvature radius, nor on the average distance of the antennae from the shower 
core position. But one of the strongest events was measured during thunderstorm activity in the vicinity of LOPES and shows the 
longest pulse length measured of 110 ns and a spectral slope of a = —3.6. 

Conclusions. We show with two different methods that frequency spectra from air shower radio emission can be reconstructed on 
event-by-event basis, with only two dozen dipole antennae simultaneously over a broad range of frequencies. According to the 
obtained spectral slopes, the maximum power is emitted below 40 MHz. Furthermore, the decrease in power to higher frequencies 
indicates a loss in coherence determined by the shower disc thickness. We conclude that a broader bandwidth, larger collecting 
■ area, and longer baselines, as will be provided by LOFAR, are necessary to further investigate the relation of the coherence, pulse 
$H ' length, and spectral slope of cosmic ray air showers. 
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1. Introduction 

Cosmic rays are particles or nuclei constantly bombard- 
ing the Earth's atmosphere with a large spread in en- 
ergy. The origin of cosmic rays with the highest energies 
is still unknown and therefore subject of inten sive research 
l|Auger Collaboration: J. Abraham et al.ll2007t ). These par- 
ticles, mostly protons, initiate an extensive shower of sec- 
ondary particles traveling with almost the speed of light 
through the air. The charged particles in the air shower 
produce electromagnetic emission relativistically beamed in 
the forward direction. This work probes the spectral de- 
pendence of the electromagnetic field produced in these 
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air showers to investigate the radio emission mechanism. 
Understanding the emission mechanism and its dependence 
on frequency and distance to the shower core will be an im- 
portant step to infer the primary particle species of the 
cosmic ray. 

The radio emission of ch arged partic l es in air showers 
was fir st observed in 1964 bvlJelley et all (|1965D at 44 MHz 
and bv lAllan fc Joneslll966l at 60 MHz. Those early studies 
found strongly pulsed radio emission from air showers de- 
tected with radio antennae triggered by particle detectors. 
A dependence of the radio emission's polarization on the 
geomagnetic fiel d favored a geom agnetic emission m echa- 
nism (see, e.g.. lAllan et al.lll967l . 1969: ISunl fl975h . The 
first quantitative radio frequency spectra were obtained 
from a few simultaneous narrowband observations below 
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100 MHz by , e.g.. ISpencerl (|l969f ) . I Allan et all ljl970T I. and 
iPrahl (|l97lf ). In this earlier work, spectra decreasing in 
frequency were obtained with spectral indices a between 
—2 and —1. However, radio measurements were abandoned 
in the late 1970s because of strong radio frequency inter- 
ference in the signal, difficulties with the interpretation of 
the measurements, and because of the success of alternative 
observing techn iques. An exc ellent review of the early work 
can be found in I Allan! l|l97lh . 

Interest in the radio technique was rev ived by the 
LOFA R Prototype Station (LOPES) group ijFalcke et all 
12005^ and by the Cosmic ray Detection Array with 
Logari thmic ElectroMagne tic Antennas (CODALEMA) 
group ijArdouin et al.l[2006l ) , both using digital data acqui- 
sition and powerful computers for data analysis on digitized 
radio signal for the first time. LOFAR is the Low Frequency 
Array. LOPES measured a linear dependence of the electric 
field amplitude of the radio emission in extensive cosmic- 
ray air showers on the energy of the primary particle. This 
result was an important step in proving the radio emission 
process to be coherent. 

The following two radio emission mechanisms are be- 
lieved to mainly contribute to the production of elec- 
tromagnetic radiation in air showers. First, the geomag- 
netic mechanism generates beamed synchrotron emission 
through relativistic electron-positron pairs in the Earth's 
magne tic field, which is believed to be the dominant mech- 
anism ijAllan et al.lll967l : lAllan et al.lll969h . This ansatz led 
to the description o f coherent geosynchrotro n emission at 
meter wavelengths ijFalcke fc Gorhaml I2003T ) . Second, rel- 
ativistic particles e mit Cherenkov r adiati on as a kind of 
photo nic Shockwave (jAskarvanl 11961 1965 : lKahn fe Lerchd 
Il966t l. However, the Cherenkov emission mechanism is more 
efficient in denser media than air, like water and ice. 
Additionally, radio emission in air showers was claimed to 
be dependent on the geoelectric field, and it was found 
that strong electric fields in thunde rstorm clouds can am- 
plify the measured radio amplitudes i j Mandolesi et al.lll974t 
iGurevich et al.ll2004 : iBuitink et al.ll2007T ). 

When a cosmic ray initiates an extensive air shower, 
a thin shower pancake of particles is formed growing to 
about a few hundred meters in diameter. The pancake has 
a thickness of a few meters, due to particle attenuation in 
the atmosphere and relativistic beaming at high gamma 
factors. Close to the axis of the cosmic-ray air shower, the 
particle pancake is thinnest and thus the pulse profile re- 
ceived on the axis is expected to be shortest. The radio 
pulse is expected to broaden with distance from the shower 
core, mainly caused by geometric-delay s in the arrival tim e 
of longitudinal shower evolution stages (jHuege et al.ll2007t ) . 
The spectral shape of the radio emission is mainly de- 
termined by the distribution of electrons and positrons 
in the shower pancake, as they emit synchrotron radia- 
tion. According to Monte Carlo simulations, the spectrum 
starts with coherent emission at wavelengths larger than 
the shower pancake thickness, and it falls off exponentially 
to higher frequencies, where coherence will be lost when the 
wavelength becomes com parable with the shower thickness 
[jHuege fc Falckell2005af ). 

A cosmic-ray radio spectrum is particularly interesting 
to probe the coherence of the geosynchrotron emission and 
to study the width of the radio pancake, which determines 
the length of the detected pulse. This is important for future 



instruments, such as LOFAR, which will provide even larger 
bandwidths and baselines than LOPES. 

We present electric field spectra determined from data 
recorded by LOPES with 30 antennae dipoles. The instru- 
ments LOPES and KASCADE are described in Sect.H The 
selection and data analysis of radio events are explained in 
Sect. [3] and in Sect. IH respectively. The obtained spectra 
are shown and discussed in Sect.[H 



2. Instrument 

The LOFAR prot otyp e station L OPES 
(www. lopes-proje ct . or g, iHorneffer et~afl I2004J : 
iFalcke et all l2005h has been optimized for the detec- 
tion of cosmic rays with energies around 10 17 eV. LOPES 
is placed in betwe en particle detectors of the experi- 
ment KASCADE i jwww- ik . f zk . de/KASCADE_home . htinll 
lAntoni et al.ll2003T ) in Karlsruhe. The data for this work 
were taken with LOPES operating 30 early LOFAR proto- 
type dipole antennae. KASCADE is the Karlsruhe Shower 
Core and Array Detector located at the Forschungszentrum 
Karlsruhe, Germany. It is made up of 16 x 16 particle 
detector stations placed on a 200 x 200 m 2 rectangular 
grid, with 13 m spacing, a central hadron calorimeter, 
and a muon tracking detector. A layout of LOPES and 
KASCADE on the site of the Forschungszentrum Karlsruhe 
is shown in Fig. Q] 

KASCADE provides LOPES with a cosmic-ray trigger 
as well as particle shower parameters, like muon number, 
electron numb er, direction of ori gin, and position of the 
shower center (jAntoni et alJ[2003 l. These parameters have 
been used to select cosmic-ray events for LOPES post- 
processing. 

The LOPES hardware digitizes the radio signals mea- 
sured by the dipole antenna wires with a rate of 80 
Megasamples per second and in a dynamic range of 12 bits. 
The frequency band used for cosmic-ray detection ranges 
from 43 MHz to 74 MHz, since the bandpass filter atten- 
uates the received signal significantly outside these limits. 
The data is stored on custom-made PC boards with a 1 GB 
ring buffer for each dipole. The readout time of a 0.8 ms 
cosmic-ray event takes about 1.5 s and the system is not 
sensitive to a new trigger for that period. The data taken 
with LOPES are matched with the KASCADE data and 
sent to storage units for further off-line analysis. LOPES 
was triggered with a large-event-trigger from KASCADE, 
which requires a detection in 10 out of the 16 detector clus- 
ters of KASCADE. 

2.1. Radio pulse broadening 

LOPES measures the pulse at antenna level, thus at dif- 
ferent positions and distances from the shower core at 
the same time. Radio frequency interference (RFI) is re- 
jected by down-weighting of narrow spectral lines that lie 
three sigmas above the average spectrum in two iterations. 
The signals of all antenna dipoles are combined by beam- 
forming to increase sensitivity and to further suppress RFI. 
The 30 antenna positions in the LOPES layout result in an 
average distance of the antennae from any possible shower 
core position within KASCADE of d = 75 m. 

Furthermore, the signal detected by the LOPES dipoles 
is bandpass filtered within the limits from 40 MHz to 
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Figure 2. Electronics gain factors are plotted for all 30 
LOPES antennae. The factors are for the correction of the 
received voltage measured in several calibration campaigns 
at KASCADE in the band from 43 MHz to 74 MHz. 



Fig ure 1. Layout o f LOPES inside the KASCADE ar- 
ray ijHornefferl l2006f ) . The 252 small squares indicate the 
KASCADE particle detectors. The triangles show the po- 
sitions of the LOPES antennae, as they were positioned for 
this work, and the circles highlight the first stage of LOPES 
with 10 operational antennae. The three rectangles mark 
the electronic stations that house the LOPES electronics. 
Each of these stations is collecting data from ten antennae 
elements. The top station holds the master clock module 
and a board for the reception of the KASCADE timestamp. 
Among the KASCADE det ectors are three KASCADE- 
Grande ijNavarra et al.l[2004T l detectors (squares) , which are 
spaced on a bigger grid to cover a larger area and thus 
higher primary particle energies. 



80 MHz. In principle, the filtered bandwidth of 40 MHz 
broadens an infinitely short pulse in time to 25 ns (At *~ 
1 / Av) and changes the spectral slope outside the band due 
to attenuation. The impulse response of a first LOPES an- 
tenna fitted with a Gaussian resulted in a width of 57 ns 
IjHornefferl 120061 ). Thus, additional broadening must be 
caused by the phase-characteristic of the LOPES hardware, 
but the spectral slope is not changed within the LOPES 
band. Therefore, the measured pulse lengths are upper lim- 
its and the minimum detectable width is close to 25 ns. 

2.2. Gain calibration 

The inverted V-shaped dipole antennae of LOPES have a 
direction dependent antenna gain and the data acquisition 
hardware as a whole has a frequency dependent electronics 
gain. The electronics gain was determined for each antenna 
in several dedicated calibration campaigns by measuring a 
well-defined signal from a reference antenna moun ted on a 
crane overhead LOPES (<Nehls et al.ll2007. subml Fig.EJ. 

The direction and frequency dependent antenna gain 
calibration values Gu t g v -\ for a LOPES antenna element, 
mounted on a m etal pedesta l, were calculated with a sim- 
ulation model by I Art si l|2005h . Here, <j> is the azimuth angle 



and 9 is the elevation angle. The pedestal and the geometry 
of the LOPES dipoles make the gain pattern frequency de- 
pendent. At the lower end of the LOPES band, the contour- 
lines of equal gain in the pattern are oval-shaped and the 
dipoles have their maximum sensitivity to the zenith with 
the major axis of the ellipse being perpendicular to the 
dipole, lying in East- West direction. With increasing fre- 
quency, the sensitivity maximum gets circular in shape con- 
centrating around the zenith (in the middle of the band) 
and then the direction of maximum gain weakens, splits, 
and moves down in elevation at azimuth angles perpendic- 
ular to the dipole (see Fig. [3|). 

The output of the analog-digital converter (ADC) can 
be converted to units of field strength pe r unit bandwidth, 
using both sets of gains ijHorneffe r 2006): 



= -T-\\-^, —Pant (1) 

Av y Gr^e,v)c 

= J_ / 47T^ ^ V]~ ^ 
y G^fi^c v Radc 

Here, E is the electric field from the air shower measured 
at the antenna; Av is the filtered bandwidth of the signal; 
v is the observing frequency; G^fi >v ) is the direction and 
frequency dependent gain of the antenna (Fig. [3|) ; P ant is 
the power received by the antenna; K v is the frequency 
dependent correction factor for the electronics; as plotted 
in Fig.[2]for each antenna; Vadc is the voltage digitized by 
the ADC; Radc is the input impedance of the ADC; Ho is 
the electromagnetic permeability; and c the speed of light. 

For earlier LOPES results, the directional gain was av- 
eraged in frequency and applied to the beam-formed time- 
series. For this work, it was applied in the frequency domain 
to each individual frequency bin (see Sect. 14. ip . 
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3. Event selection 

For this analysis, we studied events of LOPES using the 
30 antennae setup from 2005-11-16 to 2006-11-30. In this 
period, LOPES was recording in a single polarization con- 
figuration, all dipoles were parallel and most sensitive to 
radiation from the East- West direction. About one million 
events were triggered and recorded within this period. 

Events were pre-selected on parameters provided by 
KASCADE and by estimating the electric field strength: 



e cs t = (53 ±4.9) 



m MHz 



x ((1.09 ±0.017) -cos?9) 
/ -Rsa 



(2) 



x exp 



V(221 ±62)m 



10 6 / 



(0.99±0.04) 



This formula is obtained from fits to all LOPES 
events detected with the standard analysis software 
ijHorneffer et al.ll2007l ). Here, is the muon number at 
distances between 40 m and 200 m, which is provided by 
KASCADE, ■d is the geomagnetic angle (calculated using 
the shower direction reconstructed by KASCADE), and 
Rsa is the average distance of the antennae from the shower 
axis. 

We limited the distance of the shower core from the 
KASCADE center to 91 m to make sure that the shower 
core is inside the KASCADE array so that the shower re- 
construction by KASCADE provides high accuracy. For 
this work, 71 radio events were pre-selected with a lower 
limit of 4 pV m _1 MHz -1 on the electric field strength ob- 
tained with Eq. [21 

For the spectral analysis, 23 events of the 71 events 
were selected with a lower limit on the signal-to-noise ratio 
(SNR) of the E-field amplitude and the E-field root-mean- 
square (RMS) of 5 for each spectral bin. For these 23 events, 
a radio spectrum was generated with the two methods de- 
scribed in Sect.lH The estimated primary energy for these 
processed radio events lies mainly in the range from 10 17 eV 
to 10 18 eV. 

4. Data reduction 

The standard analysis software is a modular software pack- 
age and it was developed for the data reduction of LOPES 
radio events. We use this software to pre-process the data 
for this work. The software allows us to view t he data at 
any s tep in the desired sequence of processing ijHornefferl 
|2006( ). After calibration, digital filtering, and beam-forming 
of the selected events, we performed the spectral analysis 
according to the following processing pipeline. 



4.1. Processing pipeline 

Two methods have been used to produce the electric field 
spectra. First, the determination of the electric field ampli- 
tude of the cosmic-ray radio pulse detected in the antennae 
cross-correlation beam (see Sect. 14. 5p was filtered in several 
adjacent frequency bins of the LOPES band from 43 MHz 
to 74 MHz. Second, a Fourier transformation was applied 
to a few tens of samples around the full-band radio peak 
detected in the simple antennae beam (see Sect. I4.4p . 

The processing structure for the production of the 
cosmic-ray electric field spectra is shown in Fig. 01 The 



digitized signal is converted from ADC counts (01) to volt- 
age (02). A Hanning window is applied to the time-series 
data (03) to reduce leakage in the successive Fast Fourier 
Transform (FFT) (04). In the frequency domain, weights 
for the frequency dependent electronics gain (06), weights 
for correcting delays introduced by the instrument hard- 
ware (07) and weights to reduce RFI (08) are applied to 
the frequency bins (05). 

The resulting calibrated antenna signals (09) (as mea- 
sured at the antenna feeds) are multiplied with another set 
of complex weights to phase them in the direction of the cos- 
mic ray (10). The direction is provided by the KASCADE 
experiment and the standard analysis software attempts to 
refine it by fitting a maximum to the cosmic-ray pulse emis- 
sion in a radio skymap (see Sect.| 



The phased antenna spectra (12) are corrected by the 
directional gain of the antennae (11) and they are split 
in 16 sub-bands of 2.5 MHz applying modified Hanning 
window functions (13) as described in Sect. 14.21 After in- 
verse Fourier transformation, the digitally-filtered and cali- 
brated antenna time-series (14) are cross-correlated (15) to 
form the cc-beam, which is described in Sect. 14.51 For some 
events it became necessary to review the antenna selection 
and the beam forming result, and to reject antennae that 
do not contribute to a coherent pulse, which is indicated by 
the recursive arrow in Fig. [H 

From the resulting cc-beam (16) the background noise is 
subtracted (see Sect. [477]) and the amplitude of the electric 
field strength of the pulse (17) is determined by averaging 
32 samples (400 ns) around the pulse to obtain the electric 
field spectrum. 

For the second method, the phased full-band antenna 
spectra (12) are inverse Fourier transformed back to time. 
The resulting filtered and calibrated antenna time-series 
are averaged to form the simple f-beam (see Sect. l4.4p . The 
cc-beam was not used here since it calculates the abso- 
lute value of the electric field in the beam, whereas for the 
FFT-method, the full phase information is needed and thus 
the f-beam was used instead (see Sect. I4.4|l . The f-beam 
is also reduced by its background noise level and a mod- 
ified Hanning window is applied. Subsequently, the pulse 
found in the f-beam is Fourier transformed back to the fre- 
quency domain to obtain the electric field spectrum. For 
the FFT, 32 samples (400 ns) around the air shower pulse 
(see Fig. [5]) were chosen to minimize a contribution to the 
signal from the particle detectors and to obtain a spec- 
tral resolution of 2.5 MHz matching the resolution of the 
first method. The values in the top left corner of the plot 
read: in the 1st line: the event number, the LOPES fre- 
quency window from 43 MHz to 74 MHz, the estimated 
primary energy of 3.3 x 10 17 eV, the total electron number 
of 10 6 - 7 , the truncated muon number of 10 5 7 , and the geo- 
magnetic angle of 71°; in the 2nd line: the electric field am- 
plitude of the cc-beam (cc) of 15.4 pV m _1 MHz _1 , the full- 
band peak SNR of 45cr and the root-mean-square (RMS) 
of 0.3 pV m MHz ; and in the 3rd line: the electric field 
amplitude of the f-beam (f) of 15 pV nT^MHz- 1 , the full- 
band peak SNR of 70cr and the RMS of 0.2 m^MHz- 1 . 
The RMS of both full-band beams was determined on 
0.2 ms of data well before and after the air shower pulse 
(see Fig. El . 
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4.2. Frequency window function 

For the selection of the frequency sub-bands, a filter func- 
tion has been applied. The simplest window is a rectangular 
window. However, the sharp edges introduce leakage. For 
less sharp edges, a modified Hanning window is used. The 
Hanning function applies only to the beginning quarter and 
to the end quarter of the window. The Hanning window is 
basically a cosine-function and it has to be applied to the 
frequency band around the center frequency increased by 
half the period of the cosine of the Hanning window to pre- 
serve the integral of the rectangular window (see Fig. [TJ- 

4.3. Core position and beam direction 

The core position and the direction used for beam-forming 
is provided by the KASCADE experiment. The accuracy of 
the core position and shower direction of KASCADE for the 
selected showers of high energy are 1 m and 0.1 °, respec - 
tively, valid up to a zenith angle of 42° (jAntoni et al.ll2003T ) . 
The KASCADE direction is optimized by finding the posi- 
tion of maximum emission in a 4D radio image produced 
with the LOPES data. The axes of the four dimensions are 
time, azimuth angle, elevation angle, and curvature radius 
of an assumed spherical wavefront. The separation from 
the KASCADE direction for the selected events increases 
slightly with zenith angle and ranges up to 2.9° with an 
average of (1.0 ±0.5)° for zenith angles smaller than 75°, 
which is con s istent with an earlier result of (0.8 ± 0.4)° in 
iFalcke et alJ (|2005f h where only events with a maximum 
zenith angle of 60° were analyzed (one sigma statistical un- 
certainties) . 

4.4. Field strength beam 

The field strength beam (f-beam) is the average of multiple 
antennae signals being shifted by geometrical delays for the 
beam direction. The delayed antennae signals are added 
and then divided by the number of antennae: 



n). 



(3) 



Here, N is the number of antennae, / is the single antenna 
field strength time-series, and r the signal arrival time delay 
for the beam direction. 



4.5. Cross-correlation beam 

The cross-correlation beam (cc-beam) is defined as the real 
part of a standard cross-correlation of multiple antennae 
signals, shifted by one set of signal delays for a certain di- 
rection. The data from each unique antenna pair are mul- 
tiplied with each other, the resulting values are averaged, 
and then the square root is taken while preserving the sign: 



S cc {t) = ± 



\ 



N v 



N-l N 

EE/'(*-^)/i(i-^) 

i — 1 j>i 



(4) 



Here, N is the number of antennae, N pairs — N(N — 
l)/2 is the number of unique pairs of antennae, / is the 
single antenna field strength time-series, and r the signal 
arrival time delay for the beam direction. The advantage of 



the cc-beam, compared to the simple field strength beam, is 
that only the coherent signal is preserved within the beam 
and RFI, which is only seen by a single antenna is more 
strongly suppressed, since auto-correlation terms are not 
taken into account. 



4.6. Power beam 

The power beam (p-beam) is the average of all antenna 
auto-correlations: 



S,(i) 



i 



i N 



(5) 



Here, TV is the number of antennae and / is the single 
antenna field strength time-series. The power beam is sen- 
sitive to the total power received by all the antennae from 
all directions, independent of the coherence of the signal. 

4. 7. Background noise 

The system noise in a LOPES beam toward the zenith 
was determined by averaging background spectra of 664 
arbitrary chosen events recorded over a period of one year 
(11/2005 - 11/2006). For these background spectra, short 
pulses in time and narrowband signals in frequency were 
attenuated by down-weighting of single samples in time 
and single bins in frequency containing power exceeding 
one sigma above the average in two iterations. The result 
in units of electric field strength is plotted in Fig. [8] and in 
units of noise temperature is plotted in Fig. [9l The power 
in the frequency bins at 52.5 MHz and 67.5 MHz of the 
field spectrum contained strong variations, which suggests 
that not all RFI could be rejected before the calculation of 
the spectrum. The spectral field strength was obtained by 
averaging the result of the FFT applied to the filtered time- 
series of each event and the noise temperature was calcu- 
lated on the average spectral power of all events. The latter 
plot also shows the galactic noise interpolated with a model 
for 45 MHz to 408 MHz: T v = 32 K x (i//408 MHz) -2 ' 5 
(jFalcke fc Gorhamll2003r) . The galactic noise is lying only 
slightly below the system noise of LOPES. Furthermore, the 
power during two periods of 48 hours, separated by half a 
year, showed a modulation with 20% of the received signal 
in accordance with the galactic plane crossing the zenith. 
Thus, a LOPES beam is sensitive to the galactic noise, but 
still dominated by the system noise. The values outside the 
band from 45 MHz to 72.5 MHz are strongly attenuated by 
the bandpass filter. 

We calculated the level of the background noise for 
the cc-beam and for the f-beam on an event-by-event 
basis on two blocks of 0.2 ms in the first and second 
half of each event. On average in time n cc = (—0.13 ± 
0.06) fiV m^MHz- 1 and n f = (0.6 ±0.2) /iV m^MHz- 1 
(one sigma statistical uncertainty). 

The parts of the event files for the noise-calculation leave 
out the event edges, which are affected by the Hanning win- 
dow. Also, they leave out the cosmic-ray pulse and possible 
emission from the KASCADE particle detectors. 

In the LOPES data, the signal of the particle detec- 
tors is trailing the radio signal of the air shower, and 
the KASCADE trigger arrives about 1.8 /is later at the 
LOPES electronics than the radio pulse. The peak of the 
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KASCADE emission arrives shortly after the real air shower 
pulse, between 1.7 fis and 1.5 p,s before the trigger arrives 
at the LOPES electronics and it is less than 1 /is wide. The 
emission from the particle detectors is received by each an- 
tenna from several directions. Therefore, the detector noise 
does not add coherently and it is attenuated in the beam- 
forming process. 

The root-mean-square (RMS) of the background of the 
cc-beam was determined for each frequency bin on the same 
parts of each event as the noise level above. The RMS of 
the f-beam of each event was determined as an average of 
50 spectra calculated on blocks of 32 samples, which were 
chosen away from the cosmic-ray pulse. 

4.8. Slope correction 

In Fig. [TUl the slope of the 664 noise spectra recorded over 
the earlier mentioned period of one year is plotted as a 
function of zenith angle. The plot shows a steepening of 
the slopes of events with higher zenith angles. The vari- 
ation of the sky noise is averaged out over the period of 
one year and it is not observed to change with zenith an- 
gle. A small contribution might come from remaining RFI 
increasing in signal strength to the horizon. However, the 
main contribution is introduced by deviation of the fre- 
quency dependent antenna gain calibration from the real 
antenna gain. This dependence affecting the slopes of the 
air shower spectra was corrected by dividing each spectra 
by a normalized correction spectrum: 



"/v a i 



(6) 



Here, e„ is the electric field amplitude for the fre- 
quency bin v, otf is the slope parameter for the measured 
spectrum obtained from a fit to the slopes in Fig. [10] of 
-0.008 x (90 - 6/°) - 1.2, and a r is the slope of the refer- 
ence spectrum in the zenith (Fig. [8]) of —1.2. 

4.9. Pulse amplitude 

The field strength amplitude of the frequency filtered cc- 
beam was determined by averaging it over a time window 
with the size of the inverse sub-band. As an example, a 
broadened pulse for a sub-band of ~ 1.3 MHz is plotted 
in Fig. [TTJ The values in the top left corner of the plot 
read: event number 13, frequency sub-band from 48.1 MHz 
to 49.4 MHz; E-field peak SNR of 27.1<r of the frequency 
bin and the offset between the full-band peak maximum 
and the sub-band peak maximum of —20.4 ns. The dot lies 
at a higher field strength, since its frequency bin lies be- 
low the center frequency of the full-band beam of 60 MHz. 
The cc-beam shows correlation fringes for the center fre- 
quency of the sub-band. For comparison, the Gauss-fit on 
the smoothed full-band peak was broadened in width by 
the total number of frequency bins of 24. The plotted pulse 
for a sub-band below the center frequency of the LOPES 
band is higher than the full-band pulse, which is consistent 
with the negative slope in the resulting spectrum. 

4.10. Uncertainty estimation 

The uncertainty of the determined electric field strength 
is affected by the direction and frequency dependent gain 
factors, phase uncertainties introduced by the electronics 



Table 1 . List of uncertainties on the determined electric 
field strengths. 



Method 


RMS 


Gain Uncert. 


Phase Uncert. 


cc-beam 


1.59 (7.1%) 


0.35 (1.56%) 


1.13 (5%) 


f-beam 


0.78 (3.8%) 


0.32 (1.56%) 


1.01 (5%) 



These values are calculated from all 23 analyzed events in units 
of (J.V m~ MHz - . The values in brackets are the averaged per- 
centages of the actual electric field values. 



and the background noise: 



Ae 
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(7) 



Here, (|f ) is the statistical uncertainty (RMS) of the 

background noise (see Sect. 14. 7j ): (■p-) is the statistical un- 
certainty of the E-field dependence on the beam-forming 

direction; and f J^J is an estimated uncertainty of 5% in- 
troduced by pha se errors remaining after delay calibration 
(|Horneffedl2"006h . 

The dependence of the E-field on the beam-forming di- 
rection is implemented in the gain calibration and a sta- 
tistical uncertainty is estimated for each direction and fre- 
quency sub-band calculating a gradient j± on the beam- 
pattern: 



k ■ 
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Here, e is the electric field strength, k is a constant factor 
(see Eq: [l}, G (a) is the direction dependent gain factor, 
and 42 is its gradient. 

The determined uncertainties are consistent within the 
uncer tainties applicable to earlier LOPES results (jHornefferl 
l2006h . For all 23 analyzed events, the average uncertainties 
are listed in Table [TJ 



5. Results 

For the electric field spectra, a frequency resolution of 
2.5 MHz (16 bins over the 40 MHz band) was chosen as a 
compromise between pulse broadening and sufficient spec- 
tral resolution to reveal phase uncertainties. The phase 
uncertainties are mainly caused by narrowband RFI and 
deviations from the frequency dependent gain calibration. 
The radio spectra from Monte Carlo simulations wer e 
parametrized with an exponential ijHuege fc F alckc 20 05bf ) . 
For comparison, an exponential e„ = if -exp(z//MHz//3) and 
a power-law form e v = K ■ v a were fitted to the data of the 
single and the average spectra obtained by both methods 
and corrected as described in Sect. 14.71 However, the xled 
of both fit types turns out to be the same and thus both 
functions fit equally well. The resulting fit values can be 
found in Sect. 15.21 and Table [2] at the end of the paper. 
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5.1. Single event spectra 

Radio spectra were calculated for 23 LOPES events. 
Characteristic values for each event are listed in Table 
El including the event number; the event date; the 
electric field strength measured in the whole LOPES 
band from the spectra e; the spectral index a; the 
exponential fit parameter /?; the width of the Gauss- 
fitted full-band pulse At; the E-field peak SNR from 
the spectr a in sigma a; the es timated primar y en- 
ergy E v (http://www-ik.fzk.de/~r alph/CREAMl .php[ 
iGlasstetter et al.ll2005l only reliable down to 48° elevation, 
no values are given below 20° elevation); the direction (az- 
imuth <j> and elevation 0) ; the angle of the shower axis with 
the geomagnetic field the average distance of the an- 
tennae to the shower core d; and the number of antennae 
used for beam-forming N. The values for e, a, f3, At, and 
a are given as the average result from the cc-beam and the 
f-beam. The error was chosen as the maximum from both 
methods. 

Additionally, six single event spectra for both meth- 
ods are plotted in Fig. [TU For each spectrum an average 
noise spectrum was calculated by FFT of 50 blocks offset 
from the cosmic-ray radio pulse, which is plotted as well. 
Furthermore, two spectra from Monte Carlo simulations are 
added to each plot for an inclined shower with 8 = 45 ° and 
for a vertical shower with 6 = 90° [jHuege et alJl2007l ). The 
showers were simulated for a proton with a primary energy 
of E p = 10 17 eV, an azimuthal viewing angle of <f> = 45° 
from the observer measured at a shower core distance of 
100 m and were scaled to have the same average as the 
single spectra for comparison. 

The values on the top of the plots in Fig. [14] read; 
in the 1st line: the event number, the measured elec- 
tric field strength from the spectra, the E-field peak SNR 
from the spectra in sig ma, the estimated primary energy 
ijGlasstetter et alj|2005h . the width of the Gauss-fitted full- 
band pulse, the direction (azimuth angle and elevation an- 
gle) and the geomagnetic field angle; in the 2nd line: the fit 
parameters a and (3 for the spectral slope. All the E-field 
values, slope values, sigma values, and width values in the 
table and the plots in the appendix are given as the average 
of the f-beam and the cc-beam method together with the 
maximum uncertainty. 

The six plotted events have the following special char- 
acteristics. The analyzed events indicate that with larger 
pulse width, lower frequencies dominate and steepen the 
spectrum (see event [16] with a width of 51 ns and a flat 
spectrum, versus event [7] with a width of 101 ns and a 
steep spectrum) . The latter event was found to be recorded 
during thunderstorm activity in the vicinity of LOPES. 
Strong electric fields in the clouds may h ave enhanced 
the r adio emission and broadened the pulse ijBuitink et al.l 
l2007t ). However, the sample is too small and no significant 
dependence was found. Furthermore, the measured elec- 
tric field amplitude increases with increasing angle of the 
shower axis with the geomagnetic field, as expected, accord- 
ing to the geomagnetic emission mechanism (compare event 
[10] with 11 nV m-^MHz- 1 and 51°, versus event [5] with 
20 //V m _1 MHz _1 and 89°). In addition, two events with 
the largest electric field amplitude are plotted in Fig. HU 
which have the smallest uncertainties (see event [20] and 
[21] with geomagnetic angles of 78° and 70°). 



5.2. Average field spectrum 

An average field spectrum was calculated to obtain an av- 
erage spectral slope for the 23 selected LOPES events (Fig. 
fl2| . For this plot, the single event spectra were normalized 
by their mean to unity to make them energy independent 
and then they were averaged. The numbers on the top in- 
dicate the number of events taken into account for the av- 
erage spectrum and the parameters from the spectral fits 
to the plotted spectra. The average exponential fit with the 
parameter /3 is plotted for the two methods. 

The exponential parameter is obtained with /3 = 
—0.017 ± 0.004 and the spectral index is obtained with 
a = —1.0 ± 0.2. These fit parameters apply to an aver- 
age zenith angle of the 23 events of 53° and an average 
distance of the antennae from the shower core position of 
76 m. 

The av erage spectrum measured confirms basic expec- 
tations of iFalcke fc Gorhaml (|2003f h but it is not consis- 
tent within the uncertainties and it is slightly steeper than 
the slope obtained from Monte Carlo simulations based 
on air showers simulated with CORSIKA (Cosmic Ray 
Simulations for KASCADE). The fit parameters from sim- 
ulations of a 45 inclined air shower are = —0.0085 
and a = —0.49; and the values for a vertic al shower are 
= -0.012 and a = -0.7 ijHuege et al]|2007l ). For a closer 
comparison, simulations on single event basis are neces- 
sary taking into account polarization characteristics of the 
LOPES antennae and those simulations have not yet been 
performed. 

For further comparison, the spectral index recently mea- 
sured by the CODALE MA group for a sin gle event spec- 
trum is a = —1.5 ± 0.2 (jArdouin et al.ll2~006l h which is con- 
sistent with the average slope we obtained for the uncor- 
rected spectra. 

5.3. Discussion of methods 

The major difference in the determination of the electric 
field strength, between the calculation of the spectra and 
the standard analysis software, is that the latter fits a 
Gaussian to the pulse. The standard software assumes the 
pulse to be Gaussian and any deviation of the real cosmic- 
ray radio pulse shape will cause deviations in the resulting 
electric field strength. Thus, the method used here for the 
spectra is more sensitive to the shape of the radio pulses. 
The electric field strengths determined from the spectra 
agree within ±30% with those obtained with the standard 
analysis (one sigma statistical uncertainty). 

The two methods for the spectra determination are in 
statistical agreement. The noise on the cc-beam is a few per- 
cent (~3 %) larger than the noise on the f-beam (see Sect. 
14. 7p . Furthermore, the f-beam method is more accurate, 
since the statistical uncertainties of all events averaged by 
weighting is only half of the value obtained for the cc-beam 
method. Additionally, the \ 2 of the fits on the spectra from 
the f-beam method are on average an order of magnitude 
better than for the cc-beam method. 

5.4. Pulse width analysis 

The resulting average radio spectrum can be used to obtain 
an average shape for the measured radio pulses by inverse 
Fourier transformation. However, the phases are modified 
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by the electronics, which at the moment, we cannot cor- 
rect for. Thus, a lower limit on the pulse width can be 
provided by applying an inverse Fourier transformation on 
the measured spectral amplitudes (see Fig. [13]). The ob- 
tained pulse was upsampled in frequency to increase the 
number of samples for a fit. The amplitude was normal- 
ized to unity, since it is not comparable with earlier E-field 
values. A Gaussian fit on the obtained pulse results in a 
full-width-half-maximum (FWHM) of 40 ns. As expected, 
this minimum is slightly smaller than the impulse response 
of the electronics of ~ 57 ns. The average pulse width of 
the analyzed events resulted in 60 ± 20 ns (one sigma statis- 
tical uncertainty). This value was calculated on 41 events, 
for which the difference between the f-beam width and the 
cc-beam width was not larger than two samples (25 ns). 
The two samples are a quality criterion. The simulated and 
measured pulse width set an upper limit of 60 ns and the 
minimum pulse width a lower limit of 40 ns for the single 
selected LOPES pulses. Therefore, the obtained spectral 
slopes provide an upper limit, possibly broadened by un- 
known phase uncertainties, which were estimated for the 
uncertainties in the spectra. 



6. Conclusion and outlook 

With a sample of 23 strong LOPES events, we measured 
the radio spectrum received from cosmic-ray air showers 
detected with up to 24 simple dipole antennae oriented in 
the East-West direction over a frequency band of 40 MHz. 
The spectra show that a direct Fourier transformation can 
be performed on the beam-formed radio pulses measured 
with LOFAR. The accuracy with which the spectral ampli- 
tudes can be obtained is limited by the instrument noise 
and phase uncertainties. Furthermore, the quality of the 
spectral slope is limited by the quality of the antenna gain 
model, which was simulated and measured in several cali- 
bration campaigns. 

The average slope of the spectra obtained with LOPES 
confirms basic expectations, but it is not consistent within 
uncertainties and it is slightly steeper than the slope ob- 
tained from Monte Carlo simulations based on air showers 
simulated with CORSIKA. The simulations show for the 
LOPES band flatter spectra to higher zenith angles and 
larger distances from the shower core position, and they 
show steeper spectra for increasing the azimuth from direc- 
tion North to direction East. 

As expected, the spectral slopes of the selected sample 
of events depend on the length of the pulse, where longer 
pulses result in steeper spectra. However, taking into ac- 
count the low number statistics, the spectra do not show a 
significant dependence of the slope on the electric field am- 
plitude, the azimuth angle, the zenith angle, the curvature 
radius, nor on the average distance of the antennae from 
the shower core position. 

According to the obtained spectral slopes, the maximum 
power is emitted below 40 MHz. Furthermore, the decrease 
in power to higher frequencies indicates a loss in coherence 
determined by the shower disc geometry and the longitudi- 
nal distribution of particles therein. In Monte Carlo simu- 
lations, the latter was found to follow a broadened gamma- 
distribution having a long tail to several meters in shower 
thickness, which confirms a s mooth loss of cohe rence to 
higher frequencies (see Fig. 7 in lHuege et al]l2007j ). 



The selected sample of events contained one event that 
occurred during thunderstorm activity in the vicinity of 
LOPES. This thunderstorm event had a large pulse width 
and a very steep spectrum, suggesting that LOPES mea- 
sured a geoelectric effect for the first time. 

For the study of the lateral dependence of the air shower 
radio spectrum, as a function of distance to the shower axis, 
a larger array with longer baselines and more sensitive an- 
tennae is needed. LOFAR will have a very dense core with 
baselines up to 500 meters, a collecting area of ~ 10000 m 2 
at 75 MHz, and a frequency range from 10 MHz to 290 MHz 
sampled at 200 MHz, and thus will be perfectly suited to 
improve on results obtained here. Most importantly, we will 
probe the coherence of the air shower radio signal to higher 
frequencies to infer on the geometry of the shower disc. 
Understanding the shower geometry will allow us to deter- 
mine the primary particle species using radio only. 

Furthermore, within the Pierre Auger collabora- 
tion, a sparse radio antenna array is under devel- 
opme nt ijvan den Berg fc The Pierre Auger Collaboration! 
12007ft . The antennae will be tested in the field of the Pierre 
Auger Observatory (PAO, [www . auger . orgj ) . The PAO mea- 
sures cosmic-ray air showers with particle detectors in an 
area of 3000 square kilometers with distances between the 
particle detector stations of 1.5 km. Therefore, the PAO 
responds to a higher primary particle energy range than 
LOPES, which lies beyond 10 18 eV. 

A next step will be the simulation of com- 
plete air showers based on parameters provided 
by KASCADE and LOPES, including the hard- 
ware response of LOPES. Such simulations will al- 
low a close comparison of theory and experiment. 
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50 MHz 

45° 0° 315° 




AZEL Longitude 

Figure 3. Simulated gain patterns for a single LOPES 
dipole antenna, for all directions in azimuth cj> and eleva- 
tion 9, at frequencies 50 MHz (top), 60 MHz (middle) and 
70 MHz (bottom), for a bandwidth of 1 MHz. The zenith 
lies in the center of each plot and the white circles increase 
in steps of 30° in elevation down to the horizon at the plot 
edges. The antenna dipole is lying horizontally in the plots. 
The gain ranges from its minimum in black (gain = 0) to 
the maximum gain in white (gain = 5). The dark contours 
are plotted in steps of 20%. Where five contours are visible 
(middle plot), they go from 20% gain to 100% gain. 
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Figure 4. This flowchart of the cosmic-ray processing 
pipeline starts from the raw data fj [n] , where n is the index 
for the ADC samples per antenna j. The signal converted 
from ADC counts to voltage is denoted as s. The Fourier 
transformation is denoted as T and the inverse Fourier 
transformation as T^ 1 . The number of the frequency bins 
after Fourier transformation is indicated with k. Directional 
dependence of the signal and the weights is indicated by the 
vector ~ff. The outer dashed line encloses the processing 
per individual dipole and the inner dashed line encloses the 
signal calibration steps in the frequency domain. A more 
detailed description of the steps in the processing pipeline 
can be found in Sect. 14.11 
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Figure 5. Comparison of full-band cross-correlation beam 
(cc-beam, solid line), full-band field strength beam (f-beam, 
dashed), and full-band power beam (p-beam, dash dotted) 
for a characteristic strong event. Additionally, Gauss-fits to 
the cc-beam and f-beam are plotted (thin dashed line and 
thin dotted line). The vertical lines (dashed) indicate the 
time window of the f-beam for which the FFT was applied. 
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Figure 6. Comparison of the noise of the cross-correlation 
beam (cc-beam: large RMS) and the noise of the field 
strength beam (f-beam: small RMS) of a LOPES event for 
the period of 0.2 ms before and after the trigger. The peak 
slightly left of center of the time-window is the air shower 
pulse. 
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Figure 7. Modified Harming windows (dashed) and rect- 
angular windows (dotted) for three sub-bands of the fre- 
quency band from 40 MHz to 80 MHz. 
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Figure 8. System noise spectrum of the LOPES beam to 
the zenith in units of electric field strength calculated on 
664 events of 1 ms recorded over a period of one year 
(11/2005 - 11/2006). The numbers on the bottom of the 
plot are, respectively, the slope parameters a and (3 for a 
power-law function and an exponential function fitted to 
the spectrum (dashed line). 
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Figure 9. Spectrum of the LOPES beam to the zenith in 
units of noise temperature calculated on 664 events of 1 ms 
recorded over a period of one year. The numbers on the 
bottom of the plot are the slope parameters a and (3 for 
a power-law function and an exponential function fitted to 
the spectrum (dashed line). The galactic noise is plotted 
for comparison (lower dashed line). 
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Figure 10. Spectral index of 664 noise spectra as a function 
of zenith angle. The average slope is a = — 1.4 ± 0.3 and 
the fitted line is described by -0.008 x (90 - 9/°) - 1.2. 




Figure 11. Cross-correlation beam filtered in a frequency 
sub-band (thin solid line) and scaled full-band Gauss-fit 
(dashed). The vertical lines (dash-dotted) indicate the be- 
ginning and end of the part of the cc-beam that was aver- 
aged and the dot in the middle indicates the resulting value 
of 76.6 fxV m^MHz" 1 . 
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Figure 12. Comparison of average cosmic-ray electric field 
spectra obtained with 23 LOPES events by two different 
methods. Before averaging, the single spectra are normal- 
ized to unity. The frequency bin values are determined on 
the cc-beam (circles) and the Fourier transformation of the 
f-beam (squares) are fitted with one exponential, visible as 
a straight line in this log- linear plot (dashed line). 
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Figure 13. Pulse obtained by inverse Fourier transforma- 
tion of the amplitudes of the average electric field spectrum 
(solid line). The width of the Gaussian fit results in 40 ns 
(dashed line). The amplitude is normalized to unity. 



Table 2. Event parameter list for all 23 analyzed events. For more details see first paragraph of Sect. I5.ll 
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Figure 14. Single event cosmic-ray electric field spectra determined on six LOPES example events. Each plot includes 
a spectrum determined on the cc-beam (circles) for 12 sub-bands, a spectrum determined on the f-beam (squares) by 
FFT of 32 time-samples around the radio pulse, and a noise spectrum determined by FFT on 50 blocks of 32 samples 
offset from the pulse (squares in the middle of each plot). The spectra of the cosmic-ray air shower radio pulse are fitted 
with an exponential, visible as a straight line in these log- linear plots (cc-beam: dash-dotted & f-beam: dash-dot-dot). In 
addition, the simulated data of an inclined shower (dotted) and a vertical shower (dashed) are plotted. For more details 
see the second to last paragraph of Sect. 15.11 
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